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Esq., L.L.B., F.R.S., Professor of Jurisprudence in University College, London. 
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Preliminary Remarks. 

1. JLHE investigations presented in this paper consist of two parts ; the first offers 
a solution, in a certain qualified sense, of the general linear equation in finite differ- 
ences; and the second will be found to give an almost complete analysis of the 
resolution in series of the general linear differential equation with rational factors. 
The second part is deduced directly from the results of the first, although the 
subjects of which they respectively treat appear to be wholly independent of each 
other. 

With the exception of a few cases capable of solution by partial and artificial me- 
thods, there does not at present exist any mode of solving linear equations in finite 
differences of an order higher than the first ; and with reference to such equations of 
the first order, we are obliged to be content with those insufficient forms of functions 
which are intelligible only when the independent variable is an integer, and which 
may be obtained directly from the equation itself by merely giving to the indepen- 
dent variable its successive integer values. It is in this insuflScient and qualified 
sense that the solutions here given are to be taken ; and the first part of the follow- 
ing investigations may be considered as an extension of this form of solution from the 
general equation of the first order to the general equation of the nth order. 

Linear Equations in Finite Differences. 
2. A complete analytical theory of the general equation of the nth order, 

would involve its resolution into a series of equations of the first order of the form 

^x "x^'x- 1 -— ^x> 
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from which the complete solution would result, so far as known methods extend, in 
the form 

or in an expanded form, 

Mj, = \jr^ -\'"iS^a.-\ ~r ".»;-l^.»;^a!-2 v" x-i*^ x-V^ (r^ii-% \ "'\^"nXm-'r\' ••"x 

+ . . . . 

+ . . 

_1 pl(M)_i_p(»>)r<W _l_p(») pWpW _l_p(») pW p(M)/-N(«) I I .(»>)p(M)p(M) pw 

The last column contains, as is well known, the general solution of the original 
equation deprived of the terra G^ ; and the remaining columns contain the particular 
solution of the original equation. 

3. The investigation the results of which are given in this part, although it actu- 
ally succeeds in solving the original equation, will be found to contribute little or 
nothing to the analytical theory as above explained ; and this arises from the circum- 
stance, that the particular solution, instead of being produced in the separated form 
above written, is produced in an aggregated form, 

G^-fAG^_i-|-BG^_2+... ; 

and the complementary part of the general solution, instead of being produced as 
above, appears in the form of a sum of the complementary functions above written, 
the constants c', c", &c. being the same in all. But as we shall thus obtain a solution 
with an arbitrary constant, capable of solving the original equation deprived of G„ 
we have it in our power to solve the equation completely by reducing the order of 
the equation in successive steps. 

If we take the equation of the first order, 

and arrive at its solution, not properly by any analytical method, but by giving to x 
the successive values 

X, x—\, x—2 m-\-\, m, m—1, 

and eliminating all values of u^ between the first and the last, we have 

=P^(P^_.%_2+G^_0+G^ or P.P^_A_2+P^G^_,-|-G^ 

^=^"x"x-l\"x-2Ux-3'T'^x-2)~T~"x^x-l~\~^x or "x"x-l"x-2'^x-3'\~"x"x-l^x-2'rPxGx-l-\-G^ 
^(PajP*-!' •••P»j)^«i-l+(P.»;"P»!+l)Gm-|-(r^..P^+2)Gm+i+ ... +PaiPai-lG^_2-|-P^G^_l -{-G^ 
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Pp T^""T^p p T^P P 



which is the solution ordinarily given, though it is arrived at by a process somewhat 
less coarse than the above. 

By applying the same process of successive elimination to the general linear equa- 
tion (1.), and carefully observing the law by which the entrance of the factors 
P^, Q^, R^, &c. is governed, it will be found that an exactly similar solution may be 
found in the form 



u 



='-{p.-P.(2(ict:;)+0}' 



or 

where v denotes the sum of a set of distributive operations V,, Vg, Vg, . . . V„_i not of 
a strictly algebraical character, which are capable of being performed only upon fac- 
torial expressions containing consecutive values of P^, and which have the following 
significations. Vj denotes one operation of this character, signifying that the factor 
Pm-i Pm is changed into Q^ as often as it occurs, any term in which it does not occur 
disappearing, and the sum of the terms thus obtained being the result of the opera- 
tion ; so that, for example, 

v,(p._.p.)=Q., v:(p._.p.)=o, 

&c. &c. 

Again, V2 denotes another operation of a similar character, signifying that the 
factor P^-aPm-iPm is changed into R^ as often as it occurs, the result of the operation 
being as before the sum of the terms ; so that, for example, we have 

' 2\'^x—'i'^x—3'^lil—2'^x—l^3:) "^ "'x—2^x—l^x~T' '^x—i^'x—l'^ic'T' ^ic—'t'^x—S"'^} 

' 2v".<;-6 • • -"a:) ~~ '^x~3'^a:-2'^x-l"x'T' ^a:-i'^3:-2^x-l'^xl '^x-i'^x-4,":x-l"x~T' "x-5"x-4"x-3'^xi 

v:(P._a...P.)=2R._3R., &c. &C. ; 

V3 denotes the change in a similar manner of P^^aP^.jP^.jP^ into S„; 



V„_2 denotes the change in a similar manner of P„_„+2...P« into W^ ; 
and 

V„_i denotes the change in a similar manner of P„_„+i...P^ into Z^. 
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It might at first sight be supposed that this process, if successful at all, would give 
the complete complementary solution in the form 

but it will be found that in reality this introduces only one arbitrary constant. 

It now remains to place these expressions under a properly algebraic form, and to 
verify the result ; and in order to do this we must express the general term s''(P^...Pp+i) 
in terms of the factors of the original equation ; and afterwards give to p, which in 
the first instance is regarded as a constant, the successive values required for form- 
ing the several terms of the solution. 

This may be done as follows : — 

■L<6tp p — §'a;3 p p p fxi p p p p '*a;J"'P p *"« 3'IIQ p p 2v. 

Then it is easily seen that 

^ \\"x---"p + \)'^^"x'" "p-Vl\'Jp+2i(lp+Z>'--'TlxS^"x'--"p-¥l^p + \(lx-¥li 
^ ■i\"x''-^p-^-\)--^x---^p-¥\\^p-i-3'\'^p+i'T'--~T'^xS ^ "x- • • ^p+l-"p+2^x+lJ 

and generally 

iJ^P^... Pp + l) = P^. .. Pp + l { 2p + i5'^+l + Sp + s^^+l -j- 2p + 3*^+l + ... + ip+«_2M'a!+l + -ip+M- 1^(^+1 } . 

To find Vj (P^..Pp+i)3 Jt will be convenient to proceed by steps, beginning with a 
small number of terms. 
Thus 

V,(P....P._3) = (P....P._3)(y.+ ?._.+9._.), ^V^(P....P._3)=P....P._3(M.-.) 

v,{Px-Px-4)=(Px:'Px-Mq^+y.-i+q.-2+qx-s), ^v^(P,...p,_,)=p....p,_.(9.(?.-.+?._3)+^.-i9.-3) 
Vi(P^. ..P^-s) = (P... . •'Px-,)i9x+qx-i+qx-2+qx-3-hqx-4), 
^Yl{P,...P^_,)=(P,...P,_,){qXqx-2+qx-3+qx-4)+qx-i(qx-3+qx~4)+qx-2qx-4) &c. ; 

and generally, V,(P^...Pp+i) being P^...Pp+,(2;+,g'^+,), we have 

2 ' iv^ii;' ••"p+i)-—'^x'""p+i\qx'^p+iqx+iiqx-i-^p+iqx+i~T' ••'rqp+i-^p+iqx+i) 

(Similarly, it will be seen that 

^Vl{F^...Pp+r)=P^...Pp+r2%,(q,+,llll(q,,+,Tpl[q,+0)> and so on. 
In like manner we shall find 



l> ' 2V ^•••"p+lj "^•••"p+l^p+iv'il+l'^p+2^x+i) 
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2^ ^ 2\^s:---t^p+i)^^x-""p+i^p+ayx+i-^p+&v'x+i^p+ii'>"x-\-i))i &na so on, 

2 ^ m\^a:"^ p+i)^—\"ii!""p+i)^P+2m+l\Pit+l'^p+m'^it+l)i 
2 a ^ m\"x""p+l)^\"if"p+l)^p+3m+2{'^ie+l^p+2m+l\^x+l-^p+m'^x+l))) SHU SO OH, 

w^ being that term of the series qj'A, &c. which corresponds to the operation V^, as 
q^ corresponds to Vi, r^ to V^, &c. 

V,V„(P^...Pp+,)=P^...Pp+,2p+„+,+i(4+,2^;l,«>^+,), (4 corresponding to V;), 

Finally, if Nl,,p+N^,p+N"p+.. represent a series in which 
and generally 

^x.p ^ -^p+ia-liQx+l^ x-l, p) ~r -^p+Sa-lK^x+l^ x-2, p) ~T'^p+4,a-l{^a:+l^ le-s, p) '{'•■• • 
~r'^p+(n-l)a-l\^x + l'^x-n+2,p)~r-"p+na-l{^x+l^x-n+l,p) i 

then, if this series be called N^,^, we shall find that £''(P^...Pp+,) = (P^...Pp+i)N,,,p. 
It may be here observed that the number of the terms of the series Nl,p+N!^,p+.. 

cannot exceed ^{^—p)^! when a^—p is even, and ^(a?— p+1) when x—p is odd ; and 

since p has the successive values x—2, x—3, &c., it is always known whether x~p be 
odd or even. The number of terms may be less ; for if Q^ be zero, the number of 

terms would be the next whole number above ^{x—p), &c. 

That the equation 

M,=c(P,...P^+,)N,,^ 

is a complementary solution of the original equation, or in other words, a solution of 
the original equation wanting the term G^, may be directly verified ; for we have 

«.-PA-i=c(P...P^+,)(N^,p-N,_,,p) 
Now AK-i,p=0, 

^^x-l,p^qx^x-2,p~r^x^x-3,p~T'^x^x-i,p'T' • •• ~T''^x^x-n+l,ji'V^x^x-n,p} 

and generally 

AN^''+') — « N^") 4-rN^''^ -J-vN^"' -4- -|-m, N^"^ 4-5' N'"' 

'-^'■^x~l,p — 'ix^^x-2,p^^' x^^x-3,pnr*x^^x~i,pnr •• • ^U>x^^x-n+l,pn^'^x^^x-n,pi 

whence 

ANx-i^p=qi^rSx-2,p-\-^x^x-3,p~r^x^x-i,p'^ •• +^^N^-«+i,p+2;!rN,_„_p, 
and 

MDCCCL. 2 M 
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In this part of the solution, I apprehend that it will not generally be necessary to 
have regard to the lower limits, since p may have any constant value, and that value 
may be taken which is most convenient in each case; and if we make p=—cc, all 
the lower limits will have this value. 

This part of the solution would then be 

«,=cP,....P.{l+A,+B,+a+..}, 
where 

and generally each terra of the part within the brackets is formed from the preceding 
term thus : first change x into x—l, multiply by q^j^.^, and effect a summation ; then 
change (in such preceding term) a' into x—2, multiply by r^+i, and effect a summation; 
and so on until lastly we change x into x—n-\-l, multiply by z^+i, and effect a sum- 
mation ; and the sum of the parts thus obtained is the next term. 

The verification of the particular solution is easily derived from the above ; but it 
rests on the assumption that the algebraic value above given for g''(P^..Pp+i) is correct; 
to which therefore particular attention is directed. 

The equation 

«,=«"(?,.. p^+0, 

considered as a solution of the original equation wanting the term G,, implies that 

s"(P,..p^^0=P.s1P.-i--PpH-0+Q^''(P.-.--Pp+.)+RAP.-3-.P.+i)+-.+z.s"(P.-».-Pp+>)- 

Now in order to verify the particular solution, 

«*,= G. + PA-l+S''(P.P.-.)G.-.2+S''(P.P.-.P.-2)G._3+.. + S''(P...P.-;,H-0«^.-^ + -, 

it is only requisite that 

s''(P...P._^+,)=P/(P.-.-P.-^+.)+Q^^"(P-^-P-P+0+R.^''(P.-3.-P.-;>+.)+..+Z.s''(P,_„..p,_^^,); 

and this is true, for it is identical in form with the equation last above given, not- 
withstanding the occurrence of x in the lowest value of P, for this lowest value re- 
mains the same throughout the expression. 

7. If P^=0, the expressions £''(P^...P^_^) reduce themselves to those terms which 
do not contain any value of P. It would not be difficult to determine generally 
what terras these are; but probably the most convenient general method of arriving 
at the solution in an algebraical form would be to make P^ equal to a constant b, and 
to make b equal to zero in the result finally obtained. 

8. The above particular solution of the original equation is in such a form that the 
general terra of the indefinite series representing the solution is given in explicit 
terms ; but that general term may be represented in an implicit form, which perhaps 
is more convenient for practical use. 

By attending to the formation of the successive expressions e''(P^...P^), it will readily 
be seen that the series 

M A+ M.G,_, -h M,G,„, -1- . . 4- MpG,_^ -h . . 
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is a solution of the original equation, if 
M„=l, 

M,=P.Mo, 



which last is a general relation determining the coefficient of any term from the 
coefficients of the n preceding terms ; or from the coefficients of all the preceding terms 
when the number of these preceding terms is less than n. This relation is the equa- 
tion of formation, and may be regarded as universal, bearing in mind that when p 
is less than n some of the terms of this relation vanish. 

The solution of the original equation is therefore reduced to that of a similar 
equation without second member : and, by what has preceded, this solution is 

it being understood that Mo=l ; and this is the value of M^ before given. 
But the value of Mp may be found otherwise, thus : 

Make P^_p+,=Pp, Q^_j„+2=Q^, &c. &c. 

Then the solution is evidently 

M,=a''(p;..p;), 

where the operation v has the same meaning as before, except that it is applied to 
the accented letters. 
Consequently 

M,= (F,..FO{H-A,+B,+C,+ ..}, 
where 

and generally the terms are formed as stated in section 5, using p for x. 

9. I shall conclude this part of the subject with a few simple examples for the pur- 
pose of illustrating the processes here given. 

Ex. 1. Let the equation be 

That part of the solution which is independent of G^, is 

u,= c{a''+(x-l)a'-'b'+l{x-2)(x-3)a''-'b*+~{x-3){x-4)ix-5)a''-'b'+...}; 

2 M 2 
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and since in the present ease P^ and Q^ are constants, we have 

Mp=aP-]-(p-l)a''-''b^+l(p-2){p-3)aP-*b*-]-~{pS)(p~4){p-5)aP'^b'+... 
whence the particular value of u^ is a series of terms of the form 

The difference in character between the solution proposed in this paper, and that 
which would result from a perfect analysis of the general equation, may be exempli- 
fied by the present instance. 

The perfect solution of this equation is known to be 



M.« 



where a and |3 are the roots of fi—at—b^=0, and each 2 introduces a constant. 
Now if the constants be made equal, and the expression be written at length, we 
shall obtain the form derived above. 

The expansion of the two terms within the parenthesis gives 

and 

G,,,+/3G.+^U_,+..+^''+'G^+...+c/3^ 

and if the difference of these be divided by a— jS, we get 

which is the same in eflfect as the result which we have obtained. 
jBjt. 2. Let the equation be 

Then 

'^ \ I /j^ I X[X+1) ' \x{x+l\ X[X+1)J ' J 

=e|r(x+i)-«rx+^r(x-i)-^r(a^-2)+..|; 

which may be put under the form of the definite integral 

If we add a second, member G^, the equation for determining M^ is 
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whence 

from which the series for u^ can be expressed. 
Ex. 3. Let the equation be 

Then 

= c|a^+^^(a7+l)a-H2^(*+l>(«-l)(*-2)a-*+2;^(*+l)*...(*-4)a^-«+..|. 
£'x. 4. Let the equation be 

Then 

Ex. 5. Let the equation be 

Here 

„_„ N' -1 N" -^(^±11 (^ + l)(P + 2) 
y* — <*} i.^«,p — ij J-^jj.p — 2 2 ' 

N'" _ (^ + lM^-l)(^-2) (j> + 2)(j9 + l) (a^+l).r (p + 4)03 + 3)(j> + 2)(j) + l) 
«.p 2.4 2 ■ 2 "1" 2.4 ' 

^(fe) _ (a? + lMa?-l)(a?-2)(a?-3)(ar-4) _ (j) + 2)(;> + l) (a; + lMa?-l)(a?-2) 
^.i* 2.4.6 2 ■ 2.4 

, > + 4)b + 3)(j7 + 2)(j>+l) {a; + l)a;_ {p + 6){p + 5){p + 4){p+3){p + 2){p+l) 
"*" 2.4 '2 2.4.6 ' 

(„+2) (a?+l)..(a?-2OT) _ (;? + 2)(j9 + l) (a^ + l)..(a;-2OT + 2) (j> + 4)..(j? + l) (^+l)-.(a?-2w + 4) 
'■^^.p ~ 2.4..(2ot + 2) 2 ' 2.4. .2m ■" 2.4 " 2.4..(2m— 2) 

._ {p + 2m)..(p+l) (x+l)x -r- (p + 2m+2)..{p+l) 
"'"• — 2.4. .2m ■ 2 "^ 2.4..(2»J + 2) 
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The particular solution, therefore, is 

^<,=G,+P,G,_.+(l+a;)P.P,_,G,_,+2^PA-.P.-2G._3+(*H'^-2)(P...P.-3)G._, 

+ (3.r^-5^-2)(P,..P,_,)G._,+(a^^-lU'+6)(P,..P,_3)G,_8+....; 
and the part of the general solution found by this method is 

^x — C^r^. .rQJ! 1 -|- 2 ~r 2.4 "•" 2.4.6 ■"•••i* 

Ex. 6. Let the equation be 

^x'T'^l "x^x- 1 < ^i'^ x'^ x-l^x-2 I ••\ ^nV -T J. . . i:x-n+ 1 ) ^x-n ~~ ^x 5 

then the equation for determining M^ is 

Mp + oil "x-p+ 1 "^p'T(^2"x-p + l"x-p+2'^'-p-i'T'' •'T'^n\'^x-p+l- • "x-p+n)-"^p~tt^^^ " 5 

or making 

P — P' 

^ x—p + l -■■ jO? 

M^+a,p;M^_,+«,p;p;_,M^_,+..+<p;..p;_,+,)Mp_»=o ; 

the solution of which is evidently the solution of 

M^+a,M,,_,+a,Mp_,+..+a„Mp_„=0, 
multiplied by (P^..Pi). 

Let (r+a,^«-'+a,r-^+..+aJ-'=^+~g + ..., 

then M^=(P;..P;){CiAa^+»+C2Bi3^+»+..}; 

and taking the parts affected by each constant separately, it will be seen that the 
original equation reduces itself to a set of equations of the first order, 

z;^ - ofP^-i =a"~'AG^, 

^'.-^PA-I=^"-'BG,, 

so that its complete solution can be adequately represented. 

Solution of Differential Equations in Series. 

10. I now proceed to point out a method by which the processes above indicated 
may be made to give solutions of cei'tain general forms of linear differential equa- 
tions. 

In a paper on Linear Differential Equations presented by me to the Royal Society, 
and which the Society has done me the honour to publish in the Philosophical Trans- 
actions (Part L for 1848, p. 31), I have enunciated, and so far as is material to the 
present purpose, demonstrated the following theorem : — 

That if, in a linear differential expression (p{x, D)m=X and its solution u=->P{x, D)X, 
the letter x be changed into the operative symbol D and D into —x, we shall thus 
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obtain another linear differential expression <p(D, — «)m=X, the solution of which will 
be u=->p{I>,-a;)X. 

In the application of this theorem, care must be taken that the first-mentioned 
solution is so written that the operations included under the function 4' ^fe not sup- 
pressed ; and it must also be borne in mind that the expressions obtained by this in- 
terchange of symbols will not in all cases be obviously interpretable. 

For applications of this singular analytical process I beg leave to refer to the me- 
moir above cited, where it is employed for the solution, in finite terms, of extensive 
classes of linear differential equations, and equations in finite differences. So far as 
the process is legitimate, it is to be observed that it is founded on reasoning of a 
purely analytical character. It does not in any manner whatever flow from the 
calculus of operations, or depend for its validity upon the soundness of the logical 
basis on which this calculus rests. 

Now it is a remarkable property of this mechanical interchange of symbols, that it 
instantaneously converts a linear equation in finite differences into a linear differen- 
tial equation ; so that wherever the former is soluble, the latter is soluble also, pro- 
vided the result be intelligible, a condition always satisfied when the functions era- 
ployed are rational algebraical functions. 

As an instance worthy of notice, let us take the example last above given {Ex. 6.). 
Bearing in mind that M^_„=e~"''M„ the proposed interchange gives the equation (writing 
<pa; for P,), 

«-j- a,?>(D) (j'm) + a,<p(D)^(D - 1 ) (a=%) + ..+«„?>( D) .. ^(D -«+ 1 )(£«%)= G ; 

whose solution, therefore, depends upon that of 

v-a(p(D){i'v)=a''-'AG, 
a proposition established by Mr. Boole by the methods of the Calculus of Opera- 
tions. 

I propose, therefore, now to employ this theorem of the interchange of symbols for 
the purpose of converting the forms of solution, above given, of equations in finite 
differences into the particular solutions of some general forms of differential equa- 
tions ; viz. those equations whose factors do not contain any irrational or transcen- 
dental functions of a?, or contain them only in the form of series of ascending powers 

of X. 

11. Mr. Boole, in his General Method in Analysis, has shown that expressions of 
this character may be placed in the form 

by changing the independent variable from x to its logarithm 6, and making use of 
the relation, [D being ^j, 

D(D-1)..(D-«+1)m=*«(^)V 
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Making use of this relation, we immediately convert our equation, which we assume 
to be in the form 



.d^ 



^d"- 






into 



= £««G (2.) 



a„(D..(D-M+l))«+ d„((D-l)...(D-»))(s%)+ e„((D-2)...(D-»-l))(e^M+-+ U{D-p)...(D-n-p+l)){,Peu)+.. 

+«™-i((D-l).. (D-«+l))(8»«)+J„_i((D-2).. (D-n ))(A)+....+*„_i((D-i;).. (D-n-p+2)XsP»n)+.. 

+«»-2((D-2)...(D-«+l))(e2»«)+... 4-*«-2((D-jt>).. (D-n-p+3)){,P«u)+... 

+ +.. 

or 

Now if in this equation we change D into 6 and into — D, we obtain the equation 
in finite differences (which suppose to be of the nth order), 

or 

Me+ypM9_i + ^M9_2+^fl_3+-"- = (/o^)"'G9_„=He_„ (suppOSc) ; 

the solution of which, by section 8, is of the form 
where Mo=], 



,+....=(^(^-l)..(^-«+l))-G, 






Restoring the symbols, and thereby converting H^.^ into (/o(D))-'(s"*G), which call 

H, we have 

«=MoH+M,(6«H)+M,(s=''H) + .. + M™(s™''H) + ..., 

where Mo, Mi... denote a series of operations having the following significations and 
relations; — 

M„=l, 

M,+||g]Mo=0, 



FINITE DIFFERENCES AND LINEAR DIFFERENTIAL EQUATIONS. 273 

M . /.(P-i) M+^g)M=0 



/.(D-m + 1) ,fA ^-m + 2) M-D-m + r) 

""^/olD-m+l) '■^^'»-'^/o(D-»«+2) ^"-"-^^^-'^^/olD-m+r) ^^'-^m-r— " . 

or, by passing s*, s^^ &c. outside the operations by the equation 

?)(D) (s-^H) =s'"^<p(D +m)H, 
M=MoH+s«M,H+s^*M2H+..+s'»''MJi+..., (3.) 

where the general law of relation is 

Now the expression H, which is the subject of all the operations, contains n arbitrary 
constants, since/o(D) is of the «th order; that is, provided a„ is not zero. 
Let (3^ |32...|3„ be the roots ofyo^=:0 ; then the complete value of H is 

(/«(D))-(a'^G)+c/'«+c/^*+..+cy»^ ; 
of which the first terra will give us the particular solution of the original equation, 
and each of the other terms a complementary solution. Bearing in mind the equa- 
tion (p(D){i''")=(pp .iP% we see at once that the first complementary solution is 

c,5^'«(Ao+A,s«+A/«+..+AX«+..), 
or 

Ci*ft(A„+A,a?-l-A,x^+ . . + A„.a?»+ . .), 

where Aq-— 1, and the law of formation of the coefficients is 

/„(/3,+m)A„+/0,+m)A™_,+/,(|3,+m)A„_2+ ..+/,(|3i+w)A„_,=0 ; 

and the remaining complementary solutions merely require the substitution of |32.../3„ 
successively for |3„ with new constants. 

The reader will not fail to perceive the peculiarity of the series when a„ is unity, 
(which value it can always have when it is not zero) ; for in that case the roots of 
/Qt=0 are the natural numbers from to n—l inclusive, so that the series begin re- 
spectively with the terms 1, x, x^, ...x"~\ 

If a„ be zero, that is, if the factor of the highest differential coefficient of u do not 
contain an absolute terra, then in order that the transformed equation may take the 
form (2.), it will be necessary to pass s* outside the operative functions in each term, 
and divide by s*. The initial functionyo^ will then be of the form 

b„{t...{t~n + \))+a^_,{t...{t-n+2)), 
and H will be 

Similarly, if in addition to a„ being zero, we have also b„ and a„_i respectively equal 
to zero, it will be necessary to pass s^* outside the operative functions, and divide by 
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^ ; and so on if other terms of the factors are wanting. Thus we can in all cases 
obtain the transformed equation in the required shape ; but in all the cases where 
coefficients thus vanish, there is the important qualification that j\t is no longer 
necessarily of the rath order ; so that H does not necessarily contain the proper 
number of arbitrary constants. The consequences of this consideration will be after- 
wards developed ; in the meantime we proceed to consider what modifications the 
series undergo where /q^ has two or more equal roots. 

12. If there should be two roots ofyo#=0 equal to |8i, one series will be deficient ; 
and it will be supplied as follows. The expression H, or (/o(D))~'(g"*G), contains in 
that case a term of the form iti^s^'* ; and it is easily seen that 

The wanting series is, therefore, 

;;:ia;ft{log 41 +A,a;+ A,«^+ . .) +(A>+A;x^+ ..)}, 



where 



A' — ^A-P 



If there should be three roots ofyo^=:0 equal tO(3i, two series will be deficient, one 
of which will be supplied as last mentioned ; and the other by the introduction of the 
series 

k^^{ (log x)\\ + A.X+ A^a-^H- . .)+2 log <A>+ A;x^+ . .)+(A>+ A'>^+ . .)}, 

where K{=~tA &c. ; for we have a term \6V'^ ; and it is easily seen that 

And generally since, where there arep+1 equal roots /3„ we have terms 
and since 

the deficient series will be supplied by the following, taken with a different constant 
for every value of /> from unity upwards : 



hX'\ (log xf{ 1 + Ai«+ A,a?^+. .) +jo(log a;)^-'(A>+ A,ar'+. .) 

+;j^(log^)^-(A>+A;V+.,)+..+ploga.(Ar''*+A^^-''x^+..)+(Al"'^+A^"'^^+..)}. 

As a matter of convenience, when the equal roots are zero, a temporary nominal 
value should be given to them for the purpose of differentiating A„ &c, 

13. Hitherto we have attended especially to the complementary solutions, or in 
other words, regarded G as zero. The operations, however, indicated in (3.) may be 
readily performed when G is a rational function of x ; which we will suppose to be 
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cleared of fractions, and so cleared to consist of a set of terms a^x'^. Then for this 
term (/o(D))-'0'^G) =«,(/o(D))-'s<''+-)<'=«,(/o(n+p))-'2(-^)<' ; 

and we have, as before, 

«=«^*(«+^>(Bo+Bi*+B2x2+..+B„*"+..), 
where Bo=:l, and the law of formation is 

/o(»* +? + w)B^ +/i(/i +J9 + m)B„_, +/2(« -\-f + m)B^_2 + . . . = ; 
and the whole solution is found by taking all the values of j). This will undergo a 
slight alteration, as the incipient term will be of the form/o(D)(g^"""'^*G), where r factors 
of (2.) vanish by reason of some of the constant coefficients a„, &c. heing zero. 

In like manner it would be easy to represent the series if G contained log x and 
its powers ; but for most other forms it would be necessary to expand G in order 
to represent the series explicitly. The solution however is theoretically complete, 
since it consists solely in the performance of operations which are known explicit 
functions of D. 

14. Before proceeding further with the main subject, I shall illustrate this process 
by a few CKamples. 

Ex. 1. Let the equation be 

Referring to (2.), we have 

/oD=D(D-l),/iD=&2(D-l)(D-2)+ai(D-l), 

/2(D)=C2(D-2)(D-3)+6i(D-2)+ao;, 

roots of/o^=0 are and 1. 

First complementary serieSj 

Ci(Ao+V+A2a;2_|_..+A„ar"'+..), 
where 

Ao=l,andm(m-l)A,+(d2(»*-l)(»*-2)+«i(»i-l))A„_i + (c2(m-2)(m-3)+&i(m-2)+ao)A„._2=0 

is the law of formation. 

Second complementary series, 

C2a?(Ao+ Aiir+ Aga^H- .+A„x"'+ . .), 
where 

Ao=l,and(m+l)mA„+(&2(»i)(?M— l)+«i?w)A„_,4-(c2(m— l)(m— 2)+&i(?w— l)+ao)A„_3=0 
is the law of formation. 
Particular solution, 

«=^^o(D)H+s''.|-l(D)H+a^''^f2(D)H+..+s'»''^^„(D)H+.., 
where 

H=(i)'G,^|^o(I>) = l,and/o(D+m)4^„(D)+/i(D+m)4(„_,(D)+/2(D+m),^_,(D)==0 

is the law of the formation of the operative functions. 

2 N 2 



276 MR. HARGREAVE ON THE RESOLUTION OF LINEAR EQUATIONS IN 

The three laws of formation, which are here written down at length, are in sub- 
stance the same ; for from the last the others are made by merely writing the roots 
successively in lieu of D in the factors. 

Ex. 2. Let the equation be 

Referring to (2.), we have in the first instance, 

((D-2)^-n^)(s^^«)-|-p(D-2)(s^^M)+(9(D-4)D+49+r)s*''tt=G, 
whence 

roots of fat =0 are n and —n. 

First complementary series, 

c,x"{Ao+A,x+A,x^+..+A^x^+..), 
where 

Ao=l, and {{n+my—n^)A^+p{n+m)A„.,+(q{n+my-\-r)A^_^==0 

is the law of formation. 

Second complementary series, 

c^ yAf)-j- AiX -r A.2'^ -j-,,-X'A^X -\-..), 

where the law is, as before, changing the sign of n. 
Particular solution, 

where 
and 

is the law of the formation of the operative functions. 
If p=0, 9=0, and r=l, so that the equation becomes 

the law is 



whence 



A — • ■"■m-2 

"^ m{2n+m)' 



u=c,x''i^l-l(\+n)-'x'+~il+n)-\2+n)-^x*^~j^(l+nr'(2+n)~\S+nr^^^ 
+c,a:-»|l-i(l-w)-V-l-^(l-«)->(2-?2)~V-^(I-n)-'(2-«)-'(3--«)-V+..J- 



Ex. 3. Let the equation be 






FINITE DIFFERENCES AND LINEAR DIFFERENTIAL EQUATIONS. 277 

Referring to (2.), we have 

/„(D)=D^; three roots equal to zero; 

/»(D)=?; 

m'A^+5^A„_„=0, law of formation. 
Hence the complementary series are, 



+ Czjlog « (^ 1 — ^3 .r"+ ^3 ^2^)3 «'" - ^3 (2m)8 (3„)3 «;'"+ . . . j 

4-3/^11 r»_ /"i -L i-^1 -i- ^^"-L/'l -4-1-1- i-U_5 «_ r'«-4- M 

+ f^sjllog a:)' (^ 1 - ^ «"+ ^ ^ «'"— ^ (-^ ^ «'»+ . . j 

+6 log X (^- ^ a;"- (^- + 2ii j^ (2iip ^'"+ (^^ + 2^ + g;^) ^^ (2^ (3^ ^•'"+ • •) 

^ ** V «^ n^ ^ V"^ ^ (2«)^ ^ n{2n)Jn^ (2nf ■* 

./i. -1_._1_. _£_..!_ .__±_>|1_? i_^3», U. 

and the particular solution is 

H being D-'Cs^-'G). 

15. The completeness of the preceding forms of solution depends, as above inti- 
mated, upon the circumstance that the function /oCD) is of an order not lower than 
the order of the original equation. It may however be of a lower order, as would 
take place in the first of the examples above given, if ag instead of being 1 were zero 
and &2 were also zero. 

Let the equation then be {Ex. A.), 

Referring to (2.), we have, in the first instance, 

(D-l)(e''«)+(c2(D-2)(D-3)+&i(D-2)+ao)(e='''«)=G, 
or 

De*+(c2(D-l)(D-2)+6i(D-l)+ao)(8<'«)=r''G; 

so that/o(D) is of the first order, and the root is zero. 
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The only complementary function to be obtained is therefore 

«=Cijl — ao^-}-2'^o(^H-«o)'^^-2:3«o(^i+«o)(l-2c2+2&i4-.ao)a;3 

a series which is divergent, but which, as will also be seen afterwards, is finite when 
the constant coefficients are connected by the formula (/>— \)pc2-\-pbi-{-a(^=0, p being 
any positive integer. 

Further, it may happen that the series obtained by the process, even when they 
afford a complete solution in respect of the number of constants, are divergent for all 
or for some values of «■. This may evidently be the case in the solution of the second 
of the examples above given ; for the law of the coefficients, as we advance in the 
series, approximates to A^=—qA^_2. 

In these cases, other solutions in series may be obtained by resolving the equation 
in finite differences in a series of terms of the form Hg+p instead of He_p. In order 
to effect this, all that is necessary is to write ^+r for 6, (r being the order of the 
equation in finite differences), and to divide by the factor of the last term instead of 
the factor of the first term ; or in other words, we must pass e''" outside the functions 
in (2.), and multiply by s"'^; so that this equation now assumes the form 

and the equation in finite differences is 

We have now to inquire for the roots offrt=0 ; the incipient term is (/.(D))-'(«"'-'>*G), 
which call Hj ; and the particular solution will be found to be 

where 

^^„(D)=l,and/(D-m)-4.^(D)+/_,(D-m)^/,_.(D)+...=0 

is the law of the formation of the operative functions. 

The substitution of the roots off J successively for D gives the law of formation of 
the complementary series. 

Taking the last example, the transformed equation now becomes 

(c,D(D-l)+&,D+a„)M+(D+l)(s-''«)=g-^<'G. 

Let the roots of ft or cj:{t— l)+^i^+a„=0, be /3, and jSj. 
The first complementary series is 
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where Ao=l, ancl/.((3, — m)A„+((3,— m+l)A^_i=0 is the law of formation ; the series 
then are 

+C2{ similar function of jSj }, 

which is, (since /.(j3i—w)=n(w+l)c2—2«j3i—«6i), 



Cj<a^' 



-^.-1+- ^\ — _ _-__&_l_^ft-s 



2c2-2^i-Si ~2c2-2i3i-6i 2.3c2-4/3i-2V 

_____jSi /3i-l /3j-2 ^__3 1 

~2c2-2^i-6i 2.3^3- 4/3i-26i 3.4c2-6/3i-36i'^ ' +••/ 



+C2{ similar function of |32 ....}. 

When (3i or ^^ i^ ^ positive integer, one of these series is terminable ; and if both 
are positive integers, the series derived from the smaller root is terminable, and the 
other gives no result, the coefficients becoming infinite. The first of the series will 
then be found to give the same result as that produced from the divergent series 
(which is evidently terminable in form in the case indicated), except that it begins at 
the other end. In this case the other complementary solution can be found in finite 
terms by reducing the order of the equation. 

The particular solution is 

H, being (/.(D))-(r^^G). 

Let us now return to example 2, the solution of which, as above found, is in some 
cases divergent. 

The transformed equation now becomes 

(5(D+2)2+r)«+i'(D+2)(s-*«) + ((D+2)2-r2)(r^''«)=s-'«'G. 
Roots oifrt~0 are — 2+a/ (-A, which call ^^ and (Sg. 
First complementary series, 

where 

A«=l, and (y(|3,-m+2r+r)A„+^(^,-m+2)A,_.+((|3,-m+2)='-r')A„_,=:0, 

is the law of formation ; and the second complementary series is the same, using the 
other root with a different constant. 

In like manner the particular solution is easily represented. 

16. We are now in a position to discuss the character of the various series obtained 
by this process with reference to their convergency or divergency, a subject of the 
highest importance to the value of the process. The investigations which follow, 
will, it is apprehended, be found to afford a complete test of the nature of the series. 
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In general, it will of course be understood that these researches, as to convergency 
and divergency, relate only to the complementary series, or in other words, to equa- 
tions deprived of their second member ; nevertheless they will to a great extent, if 
not throughout, apply to cases in which the second member is in the form of integer 
powers of x or of log x. 

1st. When in the set of functions 

/„(D), /(D), /.(D), .../_.(D), /(D), 

the function/(D) is of a higher dimension with regard to D than any other of the 
set, or is what we shall here call the dominant function, the solution of the equation 
can always be found in a convergent series of ascending powers of a;; and if in such 
a case we solve the equation in a series of descending powers of x, which we can do 
if we please, that series is certainly divergent. 

This is immediately apparent from the consideration of the law of the coefficients 

/O,+m)A„+/(f3.+m)A„_,+...+/(/3i+»i)A._.=0 ; 

which, as m increases without limit, approaches to the form 

A -_^A --^A ^A arm A --n A -b A ~ 

assuming that all the functions are only one degree lower than A^, which is the least 
favourable case for convergency. Therefore, if the largest of the terms of the right- 
hand side of this equation be In-^^m-p, we have 

and we can therefore arrive at a point in the series at which the ratio of the coeffi- 
cient of x™ to that of x"'~'' diminishes without limit. 

It will also be observed that the series introduced by two or more equal roots are 
of the same character as the original series from which they are derived ; for, A^ being 

dh. dK 
of the form (p{^^-\-m), we have ~A^'=~r^ '■> and when A^_p is of a higher order with re- 
ference to m than A„, —^ is also of a higher order than -^, and so for the other dif- 
ferential coefficients. 

We have now merely to inquire what must be the form of the original equation 
that/o(D) may be the dominant function. Referring to (2.), we see that it is neces- 
sary that the factor of the highest diiFerential coefficient of u should contain one term 
only. If this factor be 1 or x, no restriction need be imposed on the succeeding 
factors. If it be x"^, the factor of the next lower differential coefficient must not 

contain an absolute term ; and generally, if x'' be the factor of t-^5 the factor of , „_, 

must begin with a term not lower than x''~\ that of -^^^ with a term not lower than 

x^"% and so on. 
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In short, all equations where x does not enter into the factor of ^5 and all equa- 
tions of the form 

where j? is a positive integer, are soluble in convergent series of ascending powers of ^r. 
The third example above given is an instance of this form. 
2ndly. When in the set of functions 

/(D), l_m, .../(D), /„(D), 

the function /(D) is the dominant function, the solution of the equation can always 
be found in a convergent series of descending powers of x ; and if in such a case we 
solve the equation in a series of ascending powers of x, that series is certainly diver- 
gent. 

This is apparent, as before, from the consideration of the law of the coefficients, 

/,(^,-w)A„-l-/_i(|3,-m)A„_,+...+/o(3x-?w)A„_,=0, 

which, as m increases without limit, approaches to 

in the least favourable case for convergency. 

On proceeding to inquire what must be the form of the original equation, we see 
again that it is necessary that the factor of the highest differential coefficient of u 
should contain one term only. If this be x^, then the other restrictions are, that the 
factor of the next differential coefficient must stop at w^-\ that of the next at xi-^, 
and so on. 

In short, for all integer values oip the equation 

d^u , d'^~^u d"~^u 

is soluble in a convergent series of descending powers of ^. The 4th example above 
given is an instance of this form. 
Srdly. When in the set of functions 

/o(D), /(D), .../_>(D), /(D), 

the functions /(D) and /(D) are of the same dimensions, and are both dominant 
over all the other functions, the solution of the equation can be found in a series of 
ascending powers of x, which for some values of x is convergent, and for other values 
of X is divergent ; and the solution can also be found in a series of descending powers 
of X which is divergent for all values of x for which the other series is convergent, 
and convergent for all values of x for which the other series is divergent. 

For in the ascending series the law of the coefficients approaches the form, (the 

MDCCCL. 2 o 
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coefficient of the highest power of m in/o»* being 1, and that in^m being /„,) 

A — — / A 

and in the descending series the law approaches to 

/ A — l-A 

The former series is therefore convergent for all values of x numerically less than 

{Q~'^, and divergent for all values of x numerically greater than this limit; and the 
latter series is divergent for all values of x numerically greater than this quantity, 
and divergent for all values of x numerically less. 

The equations to which this rule is applicable are of the forms, (p' being less 
than p,) 

and 

The second example above given is an instance of the first of these forms. 
4thly. When in the set of functions 

/o(D), /(D),.../_.(D), /(D), 

one or more of the intermediate functions is or are of the same order as the extreme 
functions /„(D) and /.(D), or as the highest of these two when they differ in dimen- 
sions, the series obtained by the above processes will be divergent for some values of 
X, and we have not as yet any method of deriving a convergent series corresponding 
to these values ; and if one or more of the intermediate functions be of a higher 
dimension than the extreme functions, the series obtained by the above processes will 
certainly be divergent. 

These remaining cases therefore sever into two species ; first, where some of the 
intermediate functions are of the same order as the highest of the extreme functions ; 
secondly, where one or more of the intermediate functions are dominant. 

The first of these species includes equations of the two following forms : — 

(a„x^'+&„«^'^' + -.+4^^)^+(«»-.^^-'+^»-.*^'+-0£S+--+(«o*^'-"+..0«=^ 
in which there can be no function higher than/(,(D) ; and 

in which there can be no function higher than /(D). 

In these cases it will easily be seen that the law of the coefficients of the ascending 
series, as m increases without limit, approximates to 

«„A„+&„Am-i+---+4A»_p=0, 
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and the series therefore approaches without limit to a recurring series, in which the 
constants of relation are 

_^». _f2r2 _ ^ r^ . 

and the denominator of the rational fraction, to which the residue of the series ap- 
proaches, is 

**« 
In the descending series, the law of the coefficients approximates to 

and the series approaches without limit to a recurring series, in which the constants 
of relation are 

Kyi X Oyi 1 Qiyi X. ^ 

tyi 05 tyi 3u" Vji 30" 

and the denominator of the rational fraction, to which the residue of the series ap- 
proaches, is 

When/o(D) is higher than/,(D), the ascending series alone can be used ; when/,(D) 
is higher than /o(D)j the descending series alone can be used ; and when /o(D) and 
/,(D) are of the same dimension, either may be used ; and the approximations above 
referred to render it probable that these series, notwithstanding that they may be 
divergent, are the developments of continuous algebraical expressions. 

The second of the species above referred to includes all equations which are ex- 
cluded from the preceding forms ; that is, all forms which transgress hoih the restric- 
tions to which the equation in the third case is subjected. 

Of these forms, the solutions, whether obtained in ascending or in descending series, 
are always divergent ; and the divergency appears to be of an extreme and unmanage- 
able character. In this case we have an intermediate dominant function ; and the 
convergent solutions might, from considerations of analogy, be presumed to be series 
infinite in both directions, the roots of the dominant function determining the inci- 
pient terms. 

The treatment of these forms requires the solution of the equation in finite dif- 
ferences 

..-^rX+3-HQX+2+px+,+m*+Pa-i+Qa-2+Ra-3+--=g« 

not starting from either of the two extreme terms, as is done above, but from the 
term m,, so as to get a result in the form 

M,= .. -i-M_2G,+,-|-M_.G,+,-]-MoG,-f-M,G,_,-f-MA-2+ •• 

It would probably not be difficult to show that such a solution, exists ; but I have 
not found one in a form available for the purpose to which it is desired to be applied. 

2 o2 



and 
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17. Throughout the preceding investigations, the series obtained by the processes 
here displayed undergo a modification of form in the event of the expression /„f or 
frt, as the case may be, having one or more sets of imaginary roots. 

Let there be a couple of such roots of the form 05+|3^/— 1. In the ascending 
series these roots give 

where 

/„(«+m+^y^l)A„+/(«+m+^y^)A„_j+...=0 

/„(«-fm-/3x/:r])B^+/(«+m-^y=:i)B™_.+ ...=0 
are the respective lavFS of formation. 

It is apparent, therefore, that if A„ be of the form <p{«'-\-^J^\), B^ is of the form 

Making C2=Ci, and remembering that 

^p^~i _^_^-3V-i _2 cos (|3 log X) 

^p^i_x-^^~=2v^^ sin ((3 log x), 
the sum of the two series gives the double series, 

2c,a;«(cos (i3 log«){Ao+Bo+(A,+B,>+ •■ + (A»+BJ«'»+..} 

+y~l sin (|8 log j;){A„-B„+(A,-BOx+..+(A„-B„K+..}); 

which is necessarily real, since A„+B„ is purely real, and A„ — B^ is purely imaginary. 
Making now Cj= — c^, the sum of the two series gives another double series, (making 

2/•a;«(^/^:Tcos (i3 log a;) { A„-Bo+(A -B0*+ . . +(A„-B>"'+ ..} 
4-sin(i31og^){A„+B„+(A,+B>+.. + (A„+BJ.r'"+..}), 

which is likewise real. 

The descending series may be treated in a similar manner. 

18. Most of the examples to which the preceding processes are applied have been 
taken from the paper in the Philosophical Transactions for 1844, in which Mr. Boole 
developed his new General Method in Analysis, with which the subject matter of the 
present paper is closely connected, though the methods exhibited are distinct ; unless 
indeed it should prove, that the interchange of the symbol of operation and the inde- 
pendent variable, and the general relation exhibited by Mr. Boole's fundamental 
theorem of development connecting any system of linear differential equations with 
a corresponding system of equations in finite differences, are merely different repre- 
sentations of a part of some more general method or process. 

The principal difference in results, so far as concerns the solution in series of linear 
differential equations, appears to be, that in this paper the law of relation of the 
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coefficients of each series is distinct, and substantially the same in form for all ; and 
that it is not necessary to have recourse to the method of parameters in the cases of 
equal roots, or in the case of there being a term G on the right-hand side of the 
equation. In the case of imaginary roots, the laws of the series are in the first instance 
distinct from each other, and afterwards combined in couples. 

19. The investigations contained in the latter part of this paper reduce the problem 
of the integration in finite terms of the general linear differential equation with 
rational coefficients, to the finding of an algebraical expression representing the deve- 
lopment 

AQ-\-AiX-\-A2or-\- ....-{- A^oc"'-f- ..., 

where Ao=I, and the law of relation is 

/„mA^4-/imA^_i4-/3mA^_2+ ..-}-/,mA™_,=0 ; 

the functions /o/i/2--./r being known specific functions. The series to be summed 
closely resembles a recurring series ; it diflTers from it in this particular, that the law 
of relation, instead of being constant, has a uniform and simple variation as it pro- 
gresses along the series. If the rational coefficients should be themselves infinite 
series, the process still applies, the only diffisrence being that each term would be 
formed from all the preceding terms, instead of being formed from the r immediately 
preceding terms, or from all the preceding terms when the number of them is less 
than r. 

In those cases in which the equation is soluble in finite terms by known methods, 
we are enabled to assign the algebraical expression for the series ; a result which may 
be used for the discovery of generating functions. 

Thus, taking the general equation of the first order, 

we see that 

g'^o „,+«,^+.+i,^ — A„-l-Aia?-f A2.j?^-|-..-f A^a?™-!-..., 

where Ao=l5 and the general law of the series is 

a,»«A„-f-(&,(w— I)+ao)A„_,4-(c,(m— 2)-hZ>„)A^_2+...(/i(m— n)+/Co)A»-„=0. 
If we take the general equation of the second order, 

the law of the series will be 

a,m(m-I)A^-l-(62(»»-l)(»«-2)+a,(m-l))A™.,-|-(c,(m-2)(m-3)-f6,(m-2) + ««)A,„_,+ , 

+ (4(m— w)(m— «— l)-[-Ar,(m— w)-l-^„)A^_„=0 ; 

and generally, in equations of the wth order, the law of relation involves the nth 
power of m, the number of the term sought for. Thus the determination of soluble 
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cases of linear equations resolves itself into the inquiry, in what cases can a series 
whose scale is of the nth order be resolved into a number of series whose scale is of 
the first order. 

20. I have not thought it necessary here to extend these solutions in series to 
linear partial differential equations. The process by which the extension can be 
made is well known, and has no peculiar relation to the methods here developed. 



